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Abstract. We present a formulation of general nonlinear LC circuits within the framework of Birkhof- 
fian dynamical systems on manifolds. We develop a systematic procedure which allows, under rather mild 
non-degeneracy conditions, to write the governing equations for the mathematical description of the dy- 
namics of an LC circuit as a Birkhoffian differential system. In order to illustrate the advantages of this 
approach compared to known Lagrangian or Hamiltonian approaches we discuss a number of specific 
examples. In particular, the Birkhoffian approach includes networks which contain closed loops formed 
by capacitors, as well as inductor cutsets. We also extend our approach to the case of networks which 
contain independent voltage sources as well as independent current sources. Also, we derive a general 
balance law for an associated "energy function". 

Keywords. Conservative dynamical systems, Birkhoffian differential systems, Birkhoffian vector 
fields, Electrical networks, Geometric theory. 

1 Introduction 

In this paper we give a formulation of the dynamics of LC circuits within the framework of Birkhoffian 
systems [3]. Based on the constitutive relations of the involved inductors and capacitors and on Kirchhoff's 
laws, we define a configuration space and a corresponding Birkhoffian that describes the "elementary 
work" done by a set of "generalized forces". As a matter of fact, in order to cover circuits for which the 
topological assumptions usually imposed in the literature, are not satisfied, we are forced to describe a 
single circuit by a whole family of Birkhoffian systems parameterized by a finite number of real parameters. 
Relevant values of these parameters correspond to initial values for the time evolution of certain state 
variables of the circuit. The dimension of each configuration space is given by the cardinality of a selection 
of loops that cover the whole circuit. 

In order to study the dynamics of LC circuits, various Lagrangian and Hamiltonian formulations have 
been considered in the literature (see for example [2], 0], 0, M > 

In the Lagrangian approach, a central issue is the selection of suitable coordinates and corresponding 
velocities in terms of which the Lagrangian function is expressed. A specific technique for the sometimes 
difficult task of choosing the proper Lagrangian variables is presented in paper [H] . 

More often Hamiltonian formulations have been used to describe circuit equations. In |2] it is shown how 
to construct, based on the circuit topology, canonical variables and a Hamiltonian, so that the circuit 
equations attain canonical form. 

For a more general approach including also resistors, the RLC circuits, see Brayton-Moser's approach 
In [5], under the hypothesis that the currents through the inductors and the voltages across the capacitors 
determine all currents and voltages in the circuit via Kirchhoff's law, is proved the existence of the mixed 
potential function with the aid of which the system of differential equations describing the dynamics of 
such a network is written into a special form (see §4 in 5 ) . The mixed potential function is constructed 
explicitly only for the networks whose graph possesses a tree containing all the capacitor branches and 
none of the inductive branches, that is, the network does no contain any loops of capacitors or cutsets 
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of inductors, each resistor tree branch corresponds to a current-controlled resistor, each resistor co-tree 
branch corresponds to a voltage-controlled resistor (see §13 in 0). 

In [H] , the dynamics of a nonlinear LC circuit is shown to be of Hamiltonian nature with respect to a certain 
Poisson bracket which may be degenerate, that is, non-symplectic. In this formalism, the constitutive 
relations of the inductors and capacitors are used to define the Hamiltonian function in terms of capacitor 
charges and inductor fluxes, while the topological constraints of the network graph and Kirchhoff's laws 
define the Poisson bracket on the space of capacitor charge and inductor flux variables. 
But for all those formulations, a certain topological assumption on the electrical circuit appears to be 
crucial, that is, the circuit is supposed to contain neither loops of capacitors nor cutsets of inductors. 
In CO] an d 0] the Poisson bracket is replaced by the more general notion of a Dirac structure 

on a vector space, leading to implicit Hamiltonian systems. The Hamiltonian function is the total 
electromagnetic energy of the circuit and the vectorial state space is defined by the inductors' fluxes and 
capacitors' charges. The Dirac structure on the state space is obtained from Kirchhoff's laws. In this 
formalism, it is possible to include networks which do not obey the topological assumption mentioned 
before. 

In the paper at hand we will see that the restricted class of networks involving capacitor loops and 
inductor cut sets are naturally captured by the Birkhoffian approach. We are going to discuss explicit 
examples in order to demonstrate the advantages of the Birkhoffian approach in the analysis of the 
resulting systems. Another advantage of the Birkhoffian approach is the possible inclusion of dissipative 
effects caused by resistors included in a network. It is a straight-forward matter to extend the approach 
presented here to the case of RLC circuits, that is circuits containing resistors in addition to capacitors 
and inductors. However, to start it appears to be more instructive to restrict the theory to the case of 
LC circuits. The investigation of RLC circuits will be presented into another paper. 

The following parts of the paper are organized as follows. In Section 2 we recall the basics of Birkhof- 
fian systems (see jSJ) presented from the point view of differential geometry using the formalism of jets (see 
[5] ) . Birkhoffian formalism is a global formalism of the dynamics of implicit systems of second order ordi- 
nary differential equations on a manifold. In particular, we extend the approach in [Sj to non-autonomous 
systems in order to be able to treat the case of networks with independent voltage and current sources 
later on in section 4. In Section 3, our Birkhoffian formulation of the dynamic equations of a nonlinear 
LC circuit is introduced. Properties of the corresponding Birkhoffian such as its regularity and its con- 
servativeness are also discussed in this section. For a nonlinear LC electric network each Birkhoffian of 
the family is conservative. If there exists in the network some loop which contains only capacitors the 
Birkhoffian is never regular. For such electrical circuits, we present a systematic procedure to reduce 
the original configuration space to a lower dimensional one, thereby regularizing the Birkhoffian. On the 
reduced configuration space the reduced Birkhoffian will still be conservative. In case the LC circuit has 
loops which contain only linear inductors, the original configuration space can be further reduced to a 
lower dimensional one. Inductor loops can be regarded as some conservative quantities of the network. 
In Section 4 we give a Birkhoffian formulation of a nonlinear LC circuit with independent sources and 
discuss in this context the concepts of regularity and conservativeness. For instance, it turns out that 
voltage sources do not destroy conservativeness, even in the nonlinear case, while current sources might 
do so. Finally, in Section 5 we consider some specific examples. These examples are supposed to serve 
our purpose of demonstrating the power of the Birkhoffian approach. In particular, we can allow capac- 
itor loops as well as inductor cutsets, as already mentioned before. Also, we investigate the question of 
conservativeness of the underlying Birkhoffian in case of a circuit with independent current and voltage 
sources. 

Acknowledgement. We are grateful to Professor Marsden for a fruitful discussion concerning the 
topic of this paper. 

2 Birkhoffian systems 

For a smooth m-dimensional differentiable connected manifold M, we consider the tangent bundles 
(TM,7tm,M) and (TTM, tttm, TM). Let q = (q 1 , q 2 ,..., q m ) be a local coordinate system on M. This 
induces natural local coordinate systems on TM and TTM, denoted by (q,q), respectively (q, q, dq, dq). 



2 



The 2- jets manifold J 2 (M) is a 3m-dimensional submanifold of TTM defined by 

J 2 (M) = {z G TTM / Tttm(^) = tt tm (z)} (2.1) 

where T7Tm : TTM — > TM is the tangent map of 7Tm- We write ttj :— 7Ttm\j 2 (m) = Tttm\j 2 (m)- 
(J 2 (M), ttj, TM), called the 2- jet bundle (see 0), is an affine bundle modelled on the vertical vector 
bundle (V(M), tt T m\v{M), TM), V(M) = \J veTM V v (M) , where V V (M) = {z G T V TM \ (Tttm)v(z) = 
0}. In pQ, 1 1 1 1 this bundle is denoted by T 2 (M) and named second-order tangent bundle. In natural local 
coordinates, the equality in i|2.1[l yields (g, g, d<z|j2(M)) as a local coordinate system on J 2 (M). We 
set q :— dq\j2^ M y Thus, a local coordinate system q on M induces the natural local coordinate system 
(q, q, q) on J 2 (M). For further details on this affine bundle see PP, [S], 

A Birkhoffian corresponding to the configuration manifold M is a smooth 1-form ui on J 2 (M) such 
that, for any x G M, we have 

i> = (2.2) 

where i x : f3~ 1 (x) — > J 2 (M) is the embedding of the submanifold /3 _1 (x) into J 2 (M), (3 — ttm °7Tj. From 
this definition it follows that, in the natural local coordinate system (q, q, q) of J 2 (M), a Birkhoffian uj 
is given by 

m 

w = 9, W (2-3) 

3=1 

with certain functions Qj : J 2 (M) — » R. 

The pair (M, w) is said to be a Birkhoff system (see 

The differential system associated to a Birkhoffian to (see [H] ) is the set (maybe empty) D(oj), 
given by 

L>H := {z G J 2 (M) | w(z) = 0} (2.4) 

The manifold M is the space of configurations of D(uo), and -D(w) is said to have m 'degrees of freedom'. 
The Qi are the 'generalized external forces' associated to the local coordinate system (q). In the natural 
local coordinate system, D(u>) is characterized by the following implicit system of second order ODE's 

Qj(<Z: <?7 q) = for all j = 1, m (2.5) 

We conclude that the Birkhoffian formalism is a global formalism for the dynamics of implicit systems of 
second order differential equations on a manifold. 



Let us now associate a vector field to a Birkhoffian to. 
A vector field Y on the manifold TM is a smooth function Y : TM — > TTM such that 7Ttm ° F=id. 
Any vector field 1" on TM is called a second order vector field on TM if and only if Tttm(Y v ) = v for all 
v G TM. 

A cross section A of the affine bundle (J 2 (M), ttj, TM), that is, a smooth function X : TM — > J 2 (M) 
such that 7Tj o A=id, can be identified with a special vector field on TM, namely, the second order vector 
field on TM associated to X. Indeed, because (J 2 (M), ttj, TM) is a sub-bundle of (TTM,tt tm ,TM) 
as well as of (TTM,Tttm,TM), its sections can be regarded as sections of these two tangent bundles. 
Thus, using the canonical embedding i : J 2 (M) — > TTM , X can be identified with Y", that is, Y = i o A. 
In natural local coordinates a second order vector field can be represented as 

tn 

y = E 

3=1 

A Birkhoffian vector field associated to a Birkhoffian w of M (see [S]) is a smooth second order 
vector field on TM, Y = i o A, with A : TM -> J 2 (M), such that ItoA C D(w), that is 

A*w = (2.7) 

In the natural local coordinate system, a Birkhoffian vector field is given by l|2.6|l . such that Qj(q,q,q(q,q)) = 
0. 



q 4 (q, <i) 



(2.6) 
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A Birkhoffian uj is regular if and only if 



det 



'dQj 

dq l 



{q, q, q) 



(2.8 



for all (q, q, q), and for each (g, q), there exists q such that Qj(q, q, q) = 0, j = 1, ...,m. 

If a Birkhoffian of M is regular, then it satisfies the principle of determinism, that is, there exists an 

unique Birkhoffian vector field Y = i o X associated to u> such that ImX = D(uj) (see [S]). 

A Birkhoffian uj of M is called conservative if and only if there exists a smooth function : TM —* 
R such that 

(X*u>)Y = dE u) (Y) (2.9) 
for all second order vector fields Y — i o X (see |S]). 

Equation (|2.9I) is equivalent, in the natural local coordinate system, to the identity (see [3], P- 16, eq. 4) 



J^Qjiq 9' q)q 3 = 

3=1 3=1 



9gJ 



(2.10) 



is constant on TM if and only if dE LO (Y) — for all second order vector fields Y on TM (see |Sj). 
If uj is conservative and Y is a Birkhoffian vector field, then (|2.9I) becomes 



d£ w (y) = 



(2.11) 



This means that E u is constant along the trajectories of Y. 

It is also possible to introduce, in a natural manner, the notion of constrained Birkhoff system (see 
0, §4). 

Let (M, uj) be a Birkhoff system and S a smooth constant rank affine sub-bundle of the affine bundle 
7Tj : J 2 (M) — > TM. Locally, the submanifold 6 of codimension N, is described by the vanishing of N 
independent affine functions 



b v (q, q, a) = &i (?, Q)t + a"(q, q), u = 1,N 



(2.12) 



1=1 



A triple (M, uj, 6) is called constrained Birkhoff system. 

The constrained differential system associated to the constrained Birkhoff system (M , uj, &) is 
the set 

D(uj,e) = {z e e\uj(z) = 0} (2.13) 



Let us now generalize these concepts to time-dependent dynamical systems. 
For the usual formulation of Lagrangian and Hamiltonian time-dependent mechanics (see for example 
PP, §5.1, m §4.1, §4.6 ), the configuration space has the form R x M, the phase space has the form 
R x T*M, and the velocity space has the form R x TM, with some manifold M. If (t, q) is a coordinate 
system on R x M , then (£, q, q) is a coordinate system on R x TM. Thus, R x TM can be interpreted 
as a submanifold of T(R x M) given by 

i = l (2.14) 

From the physical point of view, this means that a reference frame has been chosen. This is not the 
case for relativistic mechanics. The reference system provides a splitting between the time and the 
state coordinates of a mechanical system. Within the Birkhoffian framework, we follow the usual non- 
relativistic lines. Thus, for the time-dependent system, we have in addition the equation 

i=0 (2.15) 

In view of 12. 141 . (|2.15|) . we choose in the study of time-dependent dynamical systems the extended 
bundle R x J 2 (M). 
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A time- dependent Birkhoffian is a smooth family of 1-forms ui t on J 2 (M) denned by 



Wt = ^2Qj(t,q,q,q)dq J 
3=1 



(2.16) 



where (i, g) is the natural local coordinate system on R x J 2 (M). Thus, our time-dependent Birkhof- 
fian is obtained by merely freezing t and constructing the Birkhoffian for any fixed value of t as before. 
A time- dependent second order vector field (see |llp on R x TM has the following representation in the 
natural local coordinate system 
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3=1 



?JL+?(t, q ,q)-^ 



(2.17) 



Thus, for a time-dependent system, a time- dependent Birkhoffian vector field on R x T(M) has 
the representation (|2.17|) . where Qj(t, q, q, q(t, q, q)) = 0. 

A time-dependent Birkhoffian cut is regular if and only if 



det 



dq* 



(*> ?. q) 



(2.18) 



i,j'=l,...,m 



for all (i, g, g), and for each (t, q, q), there exists q such that Qj(t, q, q, q) = 0, j = 1, to. 

A time-dependent Birkhoffian uj t is called conservative if and only if there exists a smooth family 
of functions E UJt : TM — ► R such that, everywhere, 



J2Qj(t,q q, q)q> = ^ 

3=1 3 = 1 



9E U ■ <9£^ t .. ■ 



9gJ 



(2.19) 



If tot is conservative and Yt is a time-dependent Birkhoffian vector field then, from H2.19|l . we obtain the 
generalized balance law 



dE u 
dt 



(2.20) 



along trajectories of F t . 



3 LC circuit dynamics 

A simple electrical circuit provides us with an oriented connected graph, that is, a collection of points, 
called nodes, and a set of connecting lines or arcs, called branches, such that in each branch is given a 
direction and there is at least one path between any two nodes. A path is a sequence of branches such 
that the origin of the next branch coincides with the end of the previous one. The graph will be assumed 
to be planar, that is, it can be drawn in a plane without branches crossing. For the graph theoretic 
terminology, see, for example \7\. 

Let b be the total number of branches in the graph, n be one less than the number of nodes and to be the 
cardinality of a selection of loops that cover the whole graph. Here, a loop is a path such that the first 
and last node coincide and that does not use the same branch more than once. By Euler's polyhedron 
formula, b = m + n. 

A cutset in a connected graph, is a minimal set of branches whose removal from the graph, renders the 
graph disconnected. For example, the set branches tied to a node is a cutset. 

We choose a reference node and a current direction in each Z-branch of the graph, I = 1, ...,b. We also 
consider a covering of the graph with to loops, and a current direction in each j-loop, j — 1, to. We 
assume that the associated graph has at least one loop, meaning that to > 0. 

A graph can be described by matrices: a (6n)-matrix B G 9Jlf, n (R), rank(£?) = n, called incidence matrix 
and a (677i)-matrix A G 97tt, m (R), rank(A) = m, called loop matrix. These matrices contain only 0,1, -1. 
An element of the matrix B is if a branch b is not incident with a node n, 1 if branch b enters node 
n and -1 if branch b leaves node n, respectively. An element of the matrix A is if a branch b does not 
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belong to a loop m, 1 if branch b belongs to loop m and their directions agree and -1 if branch b belongs 
to loop m and their directions oppose, respectively. For the fundamentals of electrical circuit theory, see, 
for example 

The states of the circuit have two components, the currents through the branches, denoted by I G R b , 
and the voltages across the branches, denoted by v £ R b . Using the matrices A and B, Kirchhoff's 
current law and Kirchhoff's voltage law can be expressed by the equations 

B T I = (KCL) (3.1) 

A T v = (KVL) (3.2) 

Tellegen's theorem establishes a relation between the matrices A T and B T : the kernel of the matrix B T 
is orthogonal to the kernel of the matrix j4 T (see e.g., jS] page 5). 

The next step is to introduce the branch elements in a simple electrical circuit. The branches of the 
graph associated to an LC electrical circuit, can be classified into two categories: inductor branches and 
capacitor branches. A capacitor loop will contain only capacitor branches and an inductor cutset will 
contain only inductor branches. Let k denote the number of inductor branches and p the number of 
capacitor branches, respectively. We assume that just one electrical device is associated to each branch, 
then, we have b = k + p. Thus, we can write (i a ,i a ) £ R r x R p — R b , where i a , i a are the currents 
through the inductors, the capacitors, respectively, and v = (v a ,v a ) G R fc x R p ~ R b , where v a , v a 
describe the voltage drops across the the inductors, the capacitors, respectively. 

To exemplify, let us now write the matrices B and A, for a circuit which contains four inductors, three 
capacitors and which has the following oriented connected graph 




Figure 1 



We have fc = 4, p — 3, n — 3, m — A , b — 7 . We choose the reference node to be V4 and the 
current directions as indicated in Figure 1. We cover the graph with the loops ii, 1%, I3, I4. Let V = 
(Vi,V2,Va) G R 3 be the vector of node voltage values, 1 = (i a ,i Q ) G R 4 x R 3 be the vector of branch 
current values and v = (v a , v a ) G R 4 x R 3 be the vector of branch voltage values. 

The branches in Figure 1 are labelled as follows: the first, the second, the third and the fourth branch 
are the inductor branches Li, L2, L3, L4 and the last three branches are the capacitor branches Ci, C2, 
C3. The incidence and loop matrices, B G 9Jt73(R) and A G 37t74(R), write as 
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(3.3) 



For another choice of the covering loops and of the current directions in the loops we obtain a different 
matrix A and for another choice of the reference node and of the current directions in the branches we 
obtain a different matrix B. M 

Each capacitor is supposed to be charge-controlled. For the nonlinear capacitors we assume 

v a = C a (Q a ), a=l,...,p (3.4) 



G 



where the functions C a : R — ► R\{0} are smooth and invertible, and the Q a 's denote the charges of the 
capacitors. The current through a capacitor is given by the time-derivative of the corresponding charge 

t being the time variable. 

Each inductor is supposed to be current-controlled. For the nonlinear inductors we assume 

v a = L a (i a )^, a=l,...,k (3.6) 
where L a : R — > R\{0} are smooth invertible functions. 

If the capacitors and the inductors are linear then the relations l|3.4(l and l|3.6(l become, respectively, 

v a = — , v a = L Q — (3.7) 
c a dt 

where C Q ^ and L a ^ are distinct constants. 

Taking into account l|3.4|l . (|3.5|) . i|3.6|) . the equations (|3.1|) . I|3.2[) . become 



D 1 



la 

dt 



Ca (Qq) 



(3.8) 



In the following we give a Birkhoffian formulation for the network described by the system of equations 
Using the first set of equations we are going to define a family of m-dimensional affine- 

linear configuration spaces M c C R b parameterized by a constant vector c in R". This vector is related 
to the initial values of the Q-variables at some instant of time. At this point we notice that actually 
already the initial values corresponding to the Q-variables associated to capacitors, together with those of 
m distinguished branch currents denoted by q J below, parameterize the whole solution set of the equations 
in yd. A Birkhoffian lj c of the configuration space M c arises from a linear combination of the second 
set of equations Thus, (M c ,ui c ) will be a family of Birkhoff systems that describe the LC circuit 

considered. 

We notice that the first set of equations (|3.8|l remains exactly the same for linear and nonlinear electrical 
devices. Thus, for obtaining the configuration space, it is not important whether the devices are linear or 
nonlinear. We shall sec below that the only difference is that one ends up with a nonlinear configuration 
space or rather configuration manifold when one regularizes the resulting Birkhoffian system in the case 
of nonlinear networks. 

Let H : R 6 — > R™ be a linear map that, with respect to a coordinate system (x 1 , ...,x b ) on R b , is 
given by 

x 1 

H(x\ ...,x b ) = B T | : j (3.9) 

x b 

Then, _ff _1 (c), with c a constant vector in R", is an affine-linear subspace in R^. Its dimension is 
m = b — n, because rank(B) = n. 

We define M c as 

M c :=H-\c) (3.10) 
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We denote a coordinate system on M c by q = (q 1 , ...,q m ). Then, the natural coordinate system on the 
2-jet bundle J 2 (M C ) is (q,q,q). 

Let us now represent the Birkhoffian in a specific coordinate system on M c : 
In the vector space R fc , we identify points and vectors 

U :- (3.U, 

where (Q( a ))a=i,...,fc is a coordinate system on R fc . Taking into account l|3.11[) and the fact that the 
matrix B T is a constant matrix, we integrate the first set of equations Ij3.8[l to arrive at 

B-( Q «)=C (3.12) 

with c a constant vector in R". 

Likewise consider coordinates on R b ~ R k x R p defined by 

x 1 := Q (1) , x k := Q (fc) , x k+1 := Q u x b := Q p (3.13) 

From H3.9f) . (|3.1U|I . we see that we can define coordinates on M c by solving the equations in l|3.12[l in 
terms of an appropriate set of m of the Q-variables, say q — (q , q m ). In other words, we express any 
of the x-variables as a function of q = (q 1 , q m ), namely, 



c a = ^2 Njq 3 + const, a=l, k, 



TO 

x a =^2Af*q j + const, a = k + l,...,b (3.14) 

j=i 

with certain constants N?, and Nf- Here we can think of the constants const as being initial values of 

the a;- variables at some instant of time. 

From (E2J), (|3~TT|) . (|^7T5f) and differentiating ijTTljl we get 

I = J\fq (3.15) 

with the matrix of constants M G 50Tfc m (R), for some q £ R m . 

Using Tellegen's theorem and a fundamental theorem of linear algebra, we now find a relation between 
the matrices Af and A. By a fundamental theorem of linear algebra we have 

{KeriA 7 )) 1 - = Zm(A) (3.16) 

where A e SDT hm (R), Xer(yl T ) := {a; e R b | = 0} is the kernel of A T , Jto(j4) := {x £ R fc | Ay = 
x, for somey S R m } is the image of A and denotes the orthogonal complement in R b of the respective 
vector subspace. 

For the incidence matrix B £ 3Jlfc n (R) and the loop matrix A £ 9Jtb m (R), which satisfy Kirchhoff's law 
(|3.1|l . (|3.2|) . Tellegen's theorem writes as 

Ker{B T ) = {Ker^)) 1 - (3.17) 

From the first set of equations in (|3.8|l , and by constraction of the matrix M in l|3.15|l , we have 

Ker(B T ) = Im{N) (3.18) 

Therefore, using (|3.1t)|l . I|3.17|l . I|3.18|l . we obtain Im(A) — Im(Af). Then, another application of (|3.1(i|l 

yields 

Ker{A T ) = A'er(AA T ) (3.19) 
Taking into account (|3.19|) , we see that there exists a nonsingular matrix C £ Tl mm (R) satisfing 

CA T = M T (3.20) 
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The matrix C provides a relation between the vector of the m independent loop currents and the coordi- 
nate vector q introduced on M c . 

Taking into account fy.iyp . we define the Birkhoffian uj c of M c such that the differential system \2.5)) 
is the linear combination of the second set of equations in QcJ.ty) obtained by replacing A T with the matrix 
M T . Thus, in terms of q- coordinates as chosen before, the expressions of the components Qj(q, q, q) from 
fEM are 

Qj(<l, <7, 9) = Fj(o)Q + Gj(Q), 3 = 1) -,m 



where 



m / k 



a=l 



sl = l 



i = l 



F.m = E^ ( E-^y ) £^<r = E E^i L ° (?) ? 

i=l \ a =l / 
1 J = E MfCa-k (q) 



Gj(q) = E ^" C a-* E"^ + consi 

a=fe+l \/=l 



(3.21) 
(3.22) 
(3.23) 



a=k+l 



We claim that the Birkhoffian 13. 21]) is a conservative one. 
Indeed, for our problem, the relation (|2.10() becomes 



E 



m k 



\i=l a =l / 

or (changing the indices of summation) 

/ m k \ 



E 



vi— 1 a— 1 



E 



E 



ML 



(3.24) 



dqi dq 1 



(3.25) 



Because of the special form of the terms on the left hand side of (|3.25(l . we can look for the required 
function E Uc (q, q) as a sum of a function depending only on q, and a function depending only on q. From 
the theory of total differentials, a necessary condition for the existence of such functions is the fulfilment 
of the following relations 



dq 1 dq- 
dq 1 dq 3 







for any j, I = 1, m, where 



In view of (|3~23|) . $£27$ we get: 



i—l a—l 



dG 3 (q) 
dq 1 



a=k+l 



dm) 

dq 1 



m fe 



(3.26) 



(3.27) 



(3.28) 



(3.29) 



a=l 



= 1 a=l 



where C, 



, ._ dC a (y) J, ._ dL a ( n ) 



Ji 



Therefore, the left hand side of 1)3. 26[) become 



a=k+l 



(3.30) 
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(3.31) 



We now easily see that the expressions in (|3.3U|) . (|3.31l) are zero and (|3.26|) are satisfied. Thus, we proved 
the existence of a function E UJc (g, g) such that (|3.25|) is fulfilled. 

Let us now look for the expression of this function. For linear devices, taking into account l|3.7[l . we have 



L a (q) = L a , C a - k (q) = ^=1 — ■+ <;>nsl 



(3.32) 



Got— k 

with L a , c a being real constants. Therefore, the functions Tj{q) and Gj(q) from (|3.27|) . 1)3.23(1 become 

m. k 



i—l a—1 

b rn 



\f. a \f a 

gm)--= E E-f^ + ( co ^ 



(3.34) 

Q=fc +i 4 =i Ca ~ k 

Thus, in the linear case, it is not difficult to find the function E UJc (q,q) such that 13.2511 is satisfied. This 
is 



a—1 — l a=k+l j — 1 

In order to derive such a function for nonlinear devices, we start with the equations 



<)E, 



r = = TZl Ea=l A/?W?£„ (g) g< 



, ' E "r = G 1 ( g ) = EUi^i Q ^ (?) 



^ — "ivy; — ^ Q =fc+ 
Integrating with respect to q 1 and q , respectively, we get 

A- 



^(g 1 ,...,^,? 1 ,...,^ 1 ) = ^ [ L a (q)^ a qW?dq 1 + f 1 (q 2 ,...,q m ) + 

a=l 

b r ~ 

E C a ^ k (q)Af 1 a dq 1 +g 1 (q 2 1 ...,q m ) 
a=k+l J 

fi depends only on g 2 , <f n and gi depends only on q 2 , g m . For j = 2, we have 

ff=G2(9)=EUi^^ (<?) 
and taking into account Ij3.37|) . we obtain 

K,(g 1 ,...,« m ,« 1 ,...,n = E jLa^NttfKdf + JZa^NftfKdq 2 - 



L a (q)N1N%dq l dq 2 



+ h(q\...A m ) + 



b 

E 

a=k+l 



C a - k {q)N 1 a dq L + I C a ^ k {q)N 2 a dq 2 



C'a-k {q)M 1 a N 2 a dq 1 dq 2 



+ 3 2 (g 3 ,...,9 m ) 



(3.35) 



(3.36) 



(3.37) 



(3.38) 



(3.39) 
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which can be written in the form 

k 2 2 



a—1 l—l i\<...<ii—l 
b 2 



E M) = EE E / i. H (?^ + (i-i)i. M (?l + 



EE E (-v l+1 J Jcti\«)K---Ky^- 



.dq l > 



a=k+l 1=1 »i<...<*j=l 
/ 2 (? 3 • <? 2 (? 3 '/"') 



(3.40) 



wnere o Q .- , d „i • 

Repeating this procedure for j = 3, ... m, finally in the m — th and last step, we obtain 

k m m r r 

E u (q,q) = EE E (- 1 ) m / /[£^ 1) (^)A?f + ^# 2) (^)]^..^..^ + 

a—1 i=l ix<...<ii=l v " ^ v 

6 m m 



L/ /ft /*t n n 

EE E J Jctl\q)K---K d ^--- d ^ 



(3.41) 



a=fe+l (=1 ii<...<ii = l 



Let us now discuss the question, what to do when the Birkhoffian given by Ij3.21|l is not regular in the 
sense of definition l|2.8[l . 

If there exists at least one loop in an LC circuit that contains only capacitors, then the Birkhoffian asso- 
ciated to the network is never regular. 

Indeed, for the /-loop which contains only capacitors, on the column I of the matrix A we have Af = 
for any a = 1, k. Without loss of generality, we will assume that 1 = 1, that is 



A" = 0, for any a = 1, k 

For the Birkhoffian 13.21(1 . the determinant in (|2.8I) becomes 

k 



det 



dq H 



(?, ?, ?) 



det 



E^/^a (?) 



(3.42) 



(3.43) 



- i,j'=l,...,m 



From 



et = J2? i CfK for any a = I, k, j = 1, 
Then, taking into account 13.42[1 . we have, for example, in the case m = 2 

k 

E^/^Z a (q) = C]Cf \{A\fl x (q) + (AlfL 2 (q) + ... + {A*)*L k (q) 



Then, 



det 



EA/- Q A/?L Q (q) 



.a=l 



E(^) 2 M?) 



.a=l 



1 2 


c?c 2 


r 2 r 2 









(3.44) 



(3.45) 



since the second factor obviously vanishes. In the case m — 3, we obtain 



E^/A^M?) = C 2 C 2 



j i 



E^)^ (?) 

a=l 
k 

E( A 3) 2 ia (?) 



(cftf + c?cj) 



E^2^a (?) 



(3.46) 
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Using basic calculus, the determinant of the matrix with elements l|3.46ll can be rearranged as a linear 



combination of determinants having the columns of the form 



W °3 



piipjl 

piipji 
u 2 u 3 



u 3 ^2 
u 3 °3 



respectively, with = 2 or 3 in each case. Hence, each of those determinants contain at least two 

linearly dependent columns, that is, they vanish, and this shows that the determinant is zero in the case 
to = 3 as well. Similarly, for an arbitrary to, the determinant of the matrix with the elements 



y^\-;;\-;i ■„ ■■■): = E c 'j c '. 



ii=2 



EK) 2 M«) 



,a=l 



in 

E (fic^+qfef 



2<ii<ji 



Y,AtA a h L a {q) 



(3.47) 



is zero. 



// there exists in the network ml < to loops which contain only capacitors, all the other loops containing 
at least an inductor, we can regularize the Birkhoffian ys.21\) via reduction of the configuration 
space. The reduced configuration space M c of dimension to — ml , is a linear subspace of M c or a mani- 
fold, depending on whether the capacitors are linear or nonlinear. We claim that the Birkhoffian uj c of the 
reduced configuration space M c is still a conservative Birkhoffian. Under certain conditions on the 
functions L a , a = 1, k, which characterize the inductors, the reduced Birkhoffian ui c will be a regular 
Birkhoffian. 

Without loss of generality, we can assume that there is one loop in the network that contains only 
capacitors and in the coordinate system we have chosen 

A/"f = 0, a = 1..., k (3.48) 

Thus, the Birkhoffian components <|3.21|1 . with (|3.22|) . 13.23JI . are given by, j = 2, m, 
Ql(q,q,q) = E N?C a ^ k {q) 

rn k b 

Qj(q,q,q) = £ £>£JS/?Z„ (g) 9* + E ^^-kiq) (3.49) 

i=2 a=l a=k+l 

We note that, according to l|3.48|l . q 1 does not appear in any function L a (q) and the terms L a (g)g 1 do 
not appear in any of the Birkhoffian components <32(<Z, q,q), Qm{q, 9, <?)■ 

If the capacitors in this loop are linear devices, Q\ is a linear combination of q's and we can use this 
relation to reduce the configuration space M c , to an affine-linear subspace M c of dimension m — 1. If 
the capacitors in this loop are nonlinear devices, Q\ depends nonlinearly on the q's. We define the 
(to — l)-dimcnsional manifold M c C M c by 

b 

M c = {ge M c | E KC a -k{q) = 0} (3.50) 

Q = fc+1 

By the implicit function theorem, we obtain a local coordinate system on the reduced configuration space 
M c . Taking q 1 := q 2 ,..., g™" 1 := q m , the Birkhoffian has the form Q c = Y^=i Qjdq j , 

Qj{q,q,q) = Fj(q)q + Gj(q), where (3.51) 

rn— 1 k /m—1 \ 

FMv ■= E E-^S+D^+D^ E ^+D? ( 3 - 52 ) 

i=l a=l \ 1=1 ) 
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Gj(q) := £ ^{ j+ i)Ca-k \Kf{t, -,9 m - 1 ) + E N {i+D^ + const 



(3.53) 



a=k+l 



1=1 



f:UcK 



rn — 1 



R being the unique function such that f(qo) — q$ £ R, and 



E Kc a - k Kf(s\ •••» f 1 ' 1 ) + E -^m)^ + const = 



(3.54) 



a=k+l 



1=1 



for all q — (q 1 , ...,q m ^ 1 ) g U, with U a neighborhood of q — (q^, ..^q™^ 1 ). 

We will now prove that the Birkhoffian 13.51|l is conservative. In order to do so, we will show that there 
exists a function (q, q) satisfying 



m — 1 m— 1 

E ^ = E 

3=1 3=1 



dEu .a dE^ ~j 
-rrrrq + 
dq> Qq 3 



(3.55) 



Because of the special form of the terms on the left side of (|3.55|l . we may assume that E UJ (q,q) is a 
sum of a function depending only on q, and a function depending only on q. From the theory of total 
differentials, a necessary condition for the existence of such functions is the fulfillment of the following 
relations 



qZM _ 9Fi (g) = 

dq 1 dq 1 







9G 3 (g) _ dGi(g) = q 
dq 1 dqi 



(3.56) 



for any j, I = 1, m — 1, where 



m — 1 fe 



/ m— 1 



?M ■= E E ^ + D^ + i)^ E^(W 



q 



(3.57) 



i— 1 a— r+1 



. 1=1 



We check in the same way as for the functions J-j(q) in (|3.27|l . that the first relation in Ij3.56|l is fulfilled. 
From i|3.53|l . the second relation in (|3.56|) reads as 



E W3 + i)C' a -M 

a=k+l 



7Vl f)nl +yV ('+l) 



dq 1 

m 

dq> 



,rcc 9 M ■ ire 

^r +A/ o+i) 



(3.58) 



where C' a _ k 



:= dCa dr f n) - The relation (|3~55|) reduces to 



(3.59) 



Taking into account l|3.54[) . the above relation is fulfilled, for any j, I — 1, ...,m— 1. Indeed, taking the 
derivatives with respect to g- 7 and also to q l , in the equation (|3.54(l . we obtain, respectively, 



E X?C' a _ k {q) 

a=k+l 
b 

Y, KC' a _ k {q) 

a=k+l 



Y»df(q) ra 



1 



(3.60) 



Now we multiply in l|3.6U[l the first equation with ; the second equation with — i~- and we add the 
resulting equations to obtain the equation Q3.59J1 . 
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Thus, we proved the existence of a function E w (q, q) such that l|3.55[) is fulfilled. 

For any other loop which contains only capacitors, we just repeat this procedure. Thus, we finally arrive 
at a configuration space M c of dimension m — to', where to' denotes the total number of loops of that 
type. ■ 

In case the network has loops which contain only inductors the Birkhoffian can be a regular one 
but we can further reduce the configuration space. Inductor loops can be considered as some conserved 
quantities of the network. 

If there exists in the network to" < to loops which contain only linear inductors, all the other loops 
containing at least a capacitor, we can further reduce the configuration space. The reduced configuration 
space M c of dimension m—m", is a linear subspace of M c . We claim that the Birkhoffian cD c of the reduced 
configuration space M c is a conservative Birkhoffian. Under certain conditions on the functions L a , 
a = 1, k, which characterize the inductors, the reduced Birkhoffian uj c will be a regular Birkhoffian. 
Without loss of generality, we can assume that there is one loop in the network that contains only 
inductors and in the coordinate system we have chosen 

Af? = 0, a= l...,p (3.61) 

Thus, the Birkhoffian components (|3.21|) . with (|3.22|1 . I|3.23|l . are given by, j — 2, m, 

m k 

Qi (?,«,«) = ^^AWZ^H 1 

i—1 a—1 

m k b 

Qj{qA,q) = X^XX"^' 1 ^ E ^C a ^ k (q) (3.62) 

i=l a=l a=k+l 

We note that, according to (|3.61|l . q 1 does not appear in any function C a -k(q). 

If the inductors in this loop are linear devices, Q\ is a linear combination of ifs. We can integrate this 
relation to obtain an affine-linear relation between g's (see the first example in section 5). We can use 
this relation to reduce the configuration space M c , to an affine-linear subspace M c of dimension m — 1. 
Taking q 1 :— q 2 ,..., q m ~ l := q m , one can write the Birkhoffian components of uj c and one can prove, 
using the same ideas as in the previous reduction case, the existence of the function such that this 
Birkhoffian is conservative. 

For any other loop which contains only linear inductors, we just repeat this procedure. Thus, we finally 
arrive at a configuration space M c of dimension m — m", where m" denotes the total number of loops of 
that type. ■ 

If the devices in the m" < m inductor loops are nonlinear devices, then, Q\(q,q,q),..., Q m » (g, q, q), 
are nonlinear functions depending on g's and <j's. Using these relations we can define a smooth constant 
rank affine sub-bundle 6 C of the affine bundle irj : J 2 (M C ) — ► TM C , on which we define the constrained 
Birkhoffian system (M c , w c , 6 C ). The submanifold 6 C has codimension m" . 

4 LC electric circuits with independent current/voltage sources 

Let us now consider an electric circuit containing S/ independent current sources and Sy independent 
voltage sources, in addition to k inductors and p capacitors. Then b= k +p+S/+Sy-=m + n, where 6, to, 
n have the same meaning as in section 3. We suppose that m — Sj > 0, n — Sy > 0. The branches of the 
oriented connected graph associated to this circuit are labelled as follows: L Q , a = 1, ...,k, the inductor 
branches, C a , a = 1, ...,p, the capacitor branches, 57;, i — 1, ...,Si, the current source branches, and Sy,, 
j = l, ...,Sy, the voltage source branches. 



14 



(4.1) 



Let the basic equations governing the circuit be now written in the form 

Bf( i )+Bj( i SI (t) )=0 

V dt J 

where ^ G 9Jt (m _ S/)(fc+p )(R), ^ e aK (m _ Sl)Sv (R), flf £ 9Jl (n _ Sv)(fc+p) (R), G 9K (n _ Sv)Sj (R). We 
also assume that rank(„4f) = m — Sj, (-Bf) = n — Sy- The functions i Sl (t) and v Sv (t) are given vector 
functions of time. They describe the independent current sources and independent voltage sources, re- 
spectively. The other quantities in l|4.1|l are defined as in section 3. 

In the following we give a Birkhoffian formulation for the network described by the system of equations 
using the same procedure as in section 3. That is, using the first set of equations we are 

going to define a family of (m — Si )-dimensional affine-linear configuration spaces M c C R b parameterized 
by a constant vector c in TL n ~ Sv . A Birkhoffian oj tc on the configuration space M c arises from a linear 
combination of the second set of equations \4- Ity - Thus, (M c ,uj ta ) will be a family of Birkhoff systems 
that describe the LC circuit with independent current/voltage sources considered. 

Let £ : R k+P — > K n ~ Sv be the linear map that, with respect to a coordinate system (x 1 , ...,x k+p ) on 
R fe+P , is given by 

/ x 1 



+ Bl(i Sl {t)) (4.2) 



We define 



M c :=fr\c) (4.3) 

c being a constant vector in R"~ Sv \ M c is a time-dependent affine linear subspace in R k+P . From 
rank(£>^) = n — Sy , its dimension is k + p + Sy — n = to — Sj. 



Let us figure out the relation between i a , and coordinates on M c . As in the case without sources, 



dt ' 

taking into account l|3.11f) and the fact that the matrix Bf is a constant matrix, we integrate the first 
set of equations (|4.1|1 to arrive at 

Bl ( )+l(t) = c (4.4) 



with c a constant vector in R" Sv and I{t) a primitive of Z?J ( I SJ (t) ) . 
Likewise consider coordinates in R k+P 

x 1 := Q (1) , .., x k := Q (fe) , x k+1 := Q u .., x k+p := Q p (4.5) 

We can define coordinates on M c by solving the equations (|4.4|l in terms of an appropriate set of (to — Sj) 
of the Q- variables, say q = (q 1 , ...,q m ~ Sl ). In other words, we express any of the ^-variables as a function 
of q = (q 1 , <7 m ~ S/ ), namely, 

m— Si 

x<1 = H + /°(*) + COnst > a=1 > ■-' fc ' 

m— Si 

x a = + f a (*) + const > °t = k + l, k+p (4.6) 

3=1 

with certain constants Vt" and certain functions of t, f a (t), f a (t). 

f m a \ 

The constant matrix 71 = has rank to — Si, and there exists a nonsingular 

V jij J a=l s ...,k, a=Jc+l,...,fc+p 

j = l,...,T7l-S / 

matrix C € 9Jt( m _s 7 )( m -S/)(R) such that 

CAf = <n T (4.7) 
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We define the Birkhoffian u>t c of M c such that the differential system ]2.h^l is a linear combination of the 
second set of equations in j4-l\ ), which is obtained multiplying the second set of equations in {4-ty by the 
matrix C. Taking into account \4- 7| ), in terms of q-coordinates as chosen before, the expressions of the 
components Qj(t, q, q, q), j = 1, to — Sj are 



Q3(t,q,q,q) = Fj(t,q)q + G {t,q) + V j {t) (4.8) 



where 



<m-Si J2fai 



dt 2 , 

m-S/ / k ^ \ k ^ , 2 fa( f \ 

E E (*» ?) K + E ^ (*« ?) -^r 2 ( 4 - 9 ) 

i— 1 \ a — 1 / a—1 

k-hp / m—Sj \ fe+p 

Gj{t,q)= E ^"C—fe E ^V+TW+co^ = E WjC a - k (t,q) (4.10) 

a=fc+l V j=l J a=fe+l 

B 

V;(t) = E (C^W-v-k-P-S/W (4-11) 

sy=fc+p+S/+l 



If there exist in the network ml < m — Si loops which contain only capacitors or capacitors and indepen- 
dent voltage sources, then the Birkhoffian associated to the network is never regular. 

Indeed, in this case the functions Qj corresponding to such loops depend only on q's and t. Using the 
same procedure as in section 3, the reduced configuration space M c of dimension (m — S/) — to', is a 
linear subspace of M c or a manifold, depending on whether the devices in the loops are linear or nonlinear. 

Let us finally discuss the question whether the Birkhoffian l|4.8|l is conservative or not. 

For a linear LC circuit with independent current/voltage sources we claim that the Birkhoffian \4-£fy is 
conservative. A nonlinear LC circuit with independent current /voltage sources is conservative if and only 
if it does not contain cutsets of inductors and independent current sources. 

In order to show that the Birkhoffian (|4.8|) is conservative, we are looking for a smooth function E bJt (t, q, q) 
such that the relation i|2.19[l is fulfilled. For the Birkhoffian (|4.8fl . this relation becomes 

/ m ~ Sl k ~ \ ( k d 2 f a (t) \ 

E E (*, g) <£ hp + E m i L * ft «) + Gj{t ' q) + Vj{t) * = 

\ i=l a =l / \a=l / 

f^ + f?? (4.12) 

If the inductors and the capacitors in the network are linear devices, taking into account l|3.7f) . we easily 
find the function 



k m—Si B m— Si J\faj\fct 

1 <i<r 



E. t (t,q,q) = ^E E ^aN?Nfq l q 3 + )- E E 



a—1 ij — l a—k+l — l 



d 2 f a (t\ 

E ^L a ^l + V 3 {t) + const,}*? (4.13) 



3=1 
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which satisfies l|4.12fl . 

If the inductors and the capacitors in the network are nonlinear devices, the existence of the function 
E Wt (t,q,q) which satisfies Q4.12fl . depends on the appearance of the term X)o=i 9t°£ g (t, q) d ^ s in 
(|4.12|) . For the networks which do not contain cutsets of inductors and independent current sources, this 
term does not appear at all in i|4.12[l . In this case the proof of the existence of the function E Ut is the 
same as in the case without sources. If the term ^=1 ^j L » (*> ?) is different from zero in l|4.12|) . 

E 0^ E tt^^i ______ 

then, g 9- d .j . Therefore, the Birkhoffian 14.8|l is not conservative in the sense of definition l|2.19|) . 



5 Examples 



The first example that we present is the example from the paper (0), in which we have interchanged the 
capacitor C3 and the inductor L\ to emphasis that networks which contain capacitor loops and inductor 
cutsets fit into the formalism presented in section 3. The directed connected graph associated to this 
circuit is presented in figure 1, page [fj] 

We first suppose that all devices are linear, that is, they are described by the relations i|3.7[) . Then, 
taking into account the values of the matrices A, B given by l|3.3[) . the equations (|3.8|) which govern the 
network have the form 

t I t dQ. 2 dQ 3 

12 + 13 - si oF 

h ~ 13 - 14 = 



dt ' „ dt u 







Qi Q2 I Q3 

Ci c 2 T c 3 — u 



(5.1) 



L2 



dt 



dl 3 



j dll t dl 2 1 T 013 _ n 



-Ll 



dli 

dl 



L4 



dI L 
dl 



Ci 







where C Q ^ 0, a = 1, 2, 3 and L a ^ 0, a = 1,2,3,4, are distinct constants. The relations (|3.11|) . (|3.13|) 
read as follows for this example 



dQ( a 



(a) 



dt 



a = 1,2,3,4 

„5 



(5.2) 

x' := Q (1 ), x" := Q (4 ), x u := Q 1; x' := Q 3 (5.3) 

Using the first set of equations l|5.1|l . we define the 4-dimensional affine- linear configuration space M c . 
We solve the corresponding equations (|3.12() in terms of 4 variables. In view of the notations l|5.2|l . (|5.3() . 

we obtain, for example, 

■ const 
const 



Thus, a coordinate system on M c is given by 



The matrices of constants Af 



Af ■ 



x 1 




x 3 + 


x 5 




x 4 + 


x 6 




x 2 + 


is given 


by 




q 1 := x 7 




:= x 


Nf \ 

Mr ) m 






=1,2 


3,4, o = 




3 


=1,2,3, 





1 


1 \ 


1 











1 











1 


1 





1 


-1 1 


1 





1 





/ 



x + const 



2 q 3 



x 3 ,q 4 



(5.4) 



(5.5) 



and C in (|3~"HIJl . (f^~SD|) attain the form 



C = 



/I 




V 



(5.6) 



1 / 
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Note that, if we define the Birkhofnan uj c of M c using the second set of equations (|5.1() and not a linear 
combination of them, that is, the matrix A T instead of CA T — J\f T , then, in terms of the ^-coordinates 
introduced in (|5.5|l . we obtain u> c = Ylj=i Qji.1i 9> Q)dq 3 , with 



Qi(q,q,q) 





1 


1 




c 2 


c 3 










il H 












c 2 


c 2 



<t_ 
c 2 



c 2 



Cl 



h consi 



<3s(<?, q, q) = + (li + L 3 )<f 



const 

■4 



g 1 



const 



(5.7) 



-L ig - 3 -(L 1+ L 4 )9 4 -^-^ 
Ci Ci 

The Birkhofnan (|5.7|l is not conservative. Indeed, for the Birkhofhan (|5.7|l two of the necessary conditions 
for the existence of the function E u : TM — > R such that (|2.1()(l is fulfilled, are 



dq 2 



l 2 <t - L 2 q 

(Li + L 3 )<f 



(5.8) 



Ll<f 



Because l 2 7^ 0, we see that §^r*i 7^ f^t- Therefore, there does not exist a function £^ such that (|2.1U|) 
is fulfilled. 

However, proceeding as suggested in section 3, the functions Qj(q,q,q), j — 1,2,3,4 are given by 
(|3~?T|) . (|3~22|l . and (E]23J), that is, 



i{qA,q) 



1 L 1 , n 1 q 2 tf 
1 1 q 1 h const 



Cl 



c 2 
,1 



C3 

;2 



c 2 c 2 



Cl 



Q 2 (q, q,q) — L 2 Q 1 1 h const 



c 2 c 2 



c 2 



)s(q,q,q) 
h(q,q,q) 



(Li +L 3 )<f + Li<f 



c 2 c 2 
q 1 , 9 4 



Li<f + (li + L 4 )g 4 + — + 

Cl Ci 



c 2 
const 



const 



(5.9) 



The Birkhoffian l|5.9l) is conservative. The function E UJ (q,q) is given by (|3.35|l . that is, 

Eu(q,q) 



i Ll (<f + q 4 ) 2 + ^ 2 (q 2 ) 2 + ^ 3 (<?f + \^{ff + ^-(q 4 + q 4 ) 2 



4 

£ 



(const) jq 3 



(5.10) 



Because we are in a situation where the network has one loop which contains only capacitors, the Birkhof- 



fian corresponding to (|5.9|) is not regular. Indeed, the first row of the matrix 



9Qj 
dq> 



contains 



i,i=l,2,3,4 



only zeros, therefore det -Mf = 0. 

L aq J j,j=l,2,3,4 

As we have stated in section 3, we can reduce the configuration space from dimension 4 to dimension 3. 
Using the first equation in (|5.9|) we define M c C M c by 



{q=(q\q 2 ,q 3 ,q 4 )eM c / (- + - + 

Cl 



1 

c 2 



1 

c 3 



<£_ 

c 2 



C 2 



+ — — V const = 0} 

Cl 



(5-11) 



On the reduced configuration space M c , in the coordinate system q 1 := q 2 , q 2 := q 3 , q 3 := q 4 , the 
Birkhoffian has the form u> c = X^=i Qjdq 3 ' , 

Qi(q, q, q) = L 2 Q 1 + tiq 1 + Ci<? 2 + £ 2 q 3 + const 

Q 2 {q, q) = (li + L 3 )if + Li'f + Ci? 1 + €iq 2 + € 2 q 3 + const 

Qs{q, q, 'q) = ^if + (li + u)f + £ 2 t + £ 2 f + £ 3 g 3 + const (5.12) 
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where we have introduced the notation £1 := -i- [ 1 — 7^(7^ + 7^ + 7^- 
1 I 1 I iV 1 0* 1/1 L f J_ _i_ J_ _i_ J_ 



C2C1 \ c 



Let us now see whether the Birkhofhan (|5.12|) is regular and/or conservative. We calculate 



det 







~dQj~ 








dq 


i,j=l,2,3 



L 2 

L!+L 3 Li 
L 4 L x + L 4 



(5.13) 



Thus, if l 2 [L1L4 + l 3 (li + L4)] =/= 0, then the Birkhofhan (|5.12|) is regular. 
The corresponding Birkhofhan vector held (see section 2) , is given by: 



Y = 



,1 d 











d 



dq 



L 2 



(Li +L 3 )L 4 + L1L3 



+ l 4 ) + C2L1) g 1 + (-Ci(li + l 4 ) + C2L1) g 2 + 



(-£ 2 (li + L 4 ) + £ 3 Li) q 3 + const] 



dq 2 



(li +l 4 )l 3 + l 4 l 4 



[(-£ 2 (li + l 3 ) + dLi) q 1 + (-C 2 (li + l 3 ) + dn) g 2 + 



(-£3(1-1 + l 3 ) + £ 2 Li) (j 3 + const] 



<9<7 3 



(5.14) 



Also, the Birkhofhan (|5.12l) is conservative with the function E a (q, q) given by 

Mq,q) = \^{q+q? + \^{qf + \^{q? + \^{qf + + ?f 

1 3 

C 2 («V + ft) + 2 £ 3(9 3 ) 2 + J2( const ^ 

3=1 



(5.15) 



As we have pointed out in section 3, because the network has one loop which contains only inductors, 
this is the loop I 3 , we can further reduce the dimension of the configuration space by one. The equation 
that we use for doing this is Kirchhoff's voltage law equation for this loop, that is, the sixth equation 
in 15.1(1 . For the chosen g-coordinate system l|5.5|l and after the transformation CA T , the sixth equation 
in l|5.1|l added with the hve equation in l|5.1|l and it appeared in l|5.9|l by the function Q 3 (<7, q, q) . The 
Birkhofhan formulation was presented in a coordinate free fashion. In order to have a coordinate system 
in which the sixth equation in (|5.1I) appears in the initial form, we change the g-coordinate system by 
the following relations 



t = f 
f = f 
f = f 



(5.16) 



In terms of g-coordinates on M c , the Birkhofhan cu c — Y] ■_] Qjdq 3 , where 

Qi{q, q, q) = L2<? - ^q + ^lq 1 + £2<? 3 + const 

• •• --I "2 "3 

Q2{q,q, q) = -L 2 q + (li + l 2 + h 3 )q + h x q 

Qs(q, q, q) = L is + ( L i + L 4)<? 3 + + ^q 3 + const 

Using the second equation (|5.17l) . we dehne M c C M c by 



(5.17) 



{q = (q\q 2 ,q 3 ) G MJ - L 2 f + ( L j + L 2 + L 3 )q 2 + L^ 3 + g(t) + const = 0} 



(5.18) 
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with a certain function g(t) depending on t. 

On the reduced configuration space M c , in the coordinate system q 1 := q 1 , q 2 := q 3 , the Birkhoffian has 
the form uj = Qidq 1 + Q 2 dq 2 with 



h{q,q,q) = - — — — —9 



i(q,q, q) 



Li + L 2 + L 3 

LlL 2 

( 

Li + L 2 + L3 



Li Lo "2 „i n 

- -q +£ 1 q 1 + £ 2 q 2 



Li + L 2 + L 3 



const 



(Li + L 4 )(L 2 + L 3 ) + LjL 4 - 2 ^ ^ 2 

— <7 + £ 2 g + Z 3 q + const 



Li + L 2 + L 3 



0, ...,4, the determinant tt — . T 2 , T ^ 

' ' ' (L1+L2+L3) 



Because Lj ^ 0, z 

then, the Birkhoffian (|5.19() is regular. 

Moreover, the Birkhoffian (|5.19|l is conservative, and 



hi + L 3 Li 
LiL 2 (Li +L 4 )(L 2 + L 3 ) + L1L4 



(5.19) 

^ 0, 



Eu(q,q) = 



l 2 (li+l 3 ) .i 2 



LlL 2 



2(li +l 2 + l 3 ) 



Li + L 2 + L 3 

2 



■J. -J (Li +L 4 )(L 2 + L 3 ) + L1L4 -2 2 

9 9 H — — — — (q ) + + 



2(Li + L 2 + L 3 ) 



^(g 1 ) 2 + Crf 1 ^ + k 3 (g 2 ) 2 + ^(conat) 



j q 



(5.20) 



Let us now suppose that the inductors and the capacitors in the network are nonlinear devices, their 
constitutive relations being of the form (|3.4[) . I|3.6[) . The equations (|3.8|) which govern the network have 
now the form 

dQ 1 dQ 3 _ n 

To + To - ^2. _ ^Qa - 
12 T 13 dt dt 

Ii - 13 - 14 = 







Ci(q 1 )-C 2 (q 2 ) + C 3 (q 3 ) = 

dh 
dl 



(5.21) 



i 2 (i 2 )^ + C 2 (Q 2 )-0 



L 1 (l 1 )^-L 2 (l 2 )^+L 3 (l3)^ = 



Ci(Qi 







where L a : R — ► R\{0}, C a : R — ► R\{0} are smooth invertible functions. 

As we have pointed out in section 3, the first set of equations in (|5.21|) is the same as in the linear case, 
therefore the configuration space M c is the same, too. For the coordinate system on M c given by (|5.5I) . 
the matrices A/", C have the same expressions l|5.6f) as before. Thus, in the nonlinear case, the Birkhoffian 
becomes ui c = Ylj=i QjiiiitQ)^ w ith the functions Qj given by Q3.2ip . (|3.22|l . (|3.23l) . that is, 

Ql(q,q,q) = Csiq 1 ) + diq 1 + q 4 + const) - C 2 (-q 1 + q 2 + q 3 + const) 
Q2(q,q,q) = L 2 {q 2 )q 2 + C^-q 1 + q 2 + q 3 + const) 
Q 3 (q,q,q) = (L 3 (q 3 ) + L x {q 3 + q 4 )) if + Li(q 3 + q 4 )q 4 + 

C 2 (-V +q 2 +q 3 + const) 
Q4q,q,q) = Li(q 3 + q 4 )q 3 + (L 4 (q 4 ) + Li(q 3 + q 4 )) q 4 + Citf + q 4 + const) 



(5.22) 



The Birkhoffian l|5.22|l is conservative with the function E ul (q,q) given by i|3.41|l . that is, 



Eu(q,q) 



Li(q)(q 3 + q 4 )(dq 3 + dq 4 ) + / L 2 (q 2 )q 2 dq 2 + / L 3 (q 3 )q 3 dq 3 + 



U{q 4 )q 4 dq 4 



CMidq 1 + dq 4 ) + / C^i-dq 1 + dq 2 + dq 3 ) + / C 3 (gW 



L\(q)(<? + f)dq 3 dq 4 - / / Li(q)dq 3 dq 4 + 
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C'MdqUq 4 



C' 2 {q){-dq 1 dq 2 - dq 1 dq 3 + dq 2 dq 3 



C 2 (q)dq 1 dq 2 dq 3 



The Birkhoffian 



9Q: 

w 



i,j=l,2,3,4 



(5.23) 

contains only 



is not regular, since the first row of the matrix 

zeros. But just as in the linear case, we can reduce the configuration space from dimension 4 to dimension 
3. Different from the linear case, the reduced configuration space will not be a linear subspace of M c . 
If the functions C\, C 2l C3 are such that the Jacobian matrix for the first equation in (|5.22ll has rank 
one, we define the 3-dimensional manifold M c C M c by 

M c = {q= (q\q 2 , q 3 , g 4 ) £ M c / ^(g 1 ) + d{q l + g 4 + const) - C 2 {-q 1 + q 2 + q 3 + const) = 0} (5.24) 

By the implicit function theorem, we obtain a local coordinate system on the reduced configuration space 
M c . Taking q 1 :— q 2 , q 2 :— q 3 , q 3 := g 4 , the Birkhoffian has the form ui c — Ylj=i Qjdcp , with 



Qi(q,q,q) 

Q2(q, q,q) 

Q3(q, g,g~) 



L 2 (q)q + C 2 (-f{q\q 2 , q 3 ) + g 1 + g 2 + const) 



= (L 3 (t) + L X (r + r) ) q + Li(tt + frf + 



C 2 (~f(q 1 ,q 2 , f) +?+? + const) 



Ll(q + q)q + (U(q) + L x (q + q) ) f + 
C 1 (f(q\q 2 ,q 3 ) + q 3 + const) 



(5.25) 



/ : U C R 3 — ► R 1 being an unique function such that /(go) = Qo, Qo £ R, and C 3 (f(q)) + Ci{f(q) + g 3 + 
const) -C 2 (-f(q) + q 1 +q 2 + const) = 0, Vg = (g^g 2 ,*? 3 ) E U, with [/ a neighborhood of g = (go,g 2 ,go). 
On account of Li, L 2 , L 3 , L4 : R — ► R\{0}, we have 



L2(q) 













Lstf) + Litf + f) L^+q) 
L^ + q 3 ) i 4 (g 3 ) + ii(g 2 + g 3 ) 



^0 



(5.26) 



then, the Birkhoffian (|5.25() is regular. 

Because the network has one loop which contains only inductors, let us perform a further reduction of the 
dimension of the configuration space by one, just as we have done in the linear case. In the coordinate 
system q defined in (|5.16|) . the Birkoffian Q c — Qjdq 3 , where 



Qi{q,q, q) 



L 2 (q ~q)q ~L 2 (q - q )q + C 2 (g\ g 2 , q 6 ) 



-L 2 (q 



q)q + 



Li{q + q 3 ) + L 2 (q - q) + L 3 (q) 



"2 

q + 



t f- 2 1 - 3 V- 3 
Li{q +q)q 

!,3 S "2 



3 (g,g,g) = L t {q + q )q + L x {q + q ) + L 4 (g ) g + C^g 1 , q\ q 3 ) 



(5.27) 



Using the second equation in (|5.27(l . we can define a smooth constant rank affine sub- bundle S c of the 
affine bundle ttj : J 2 (M C ) — ► TM C via 



6 C = {($,U)e J 2 (Mc)/ -L 2 (il 1 ~<f)t+ [L^ + ^ + L^ 1 -if) + L 3 (if ) 

/.,<:<'/ ■ q)t '»} 

The constraint (5 C is integrable, in the sense that we have the foliation 

- • 1 • 2 -2-3 -2 

S const {(q,q) £T(M C )/ - £ 2 (q ~q)+C 1 (q + q ) + C 3 (q )= const} 



(5.28) 



(5.29) 
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Thus, in the nonlinear case, we draw the conclusion that we can further reduce the configuration space 
only if it is possible to find from (|5.29(l new configuration coordinates q 1 , q 2 , that is, when the constraint 
()5.29(l is holonomic. 



In order to underline that, depending on the topology of the networks with independent sources, the 
associated Birkhoffian is conservative or not, we consider the circuit shown in figure 2 below. This circuit 
contains a loop formed by capacitors and independent voltage sources and a cutset formed by inductors 
and independent current sources. We shall see that the Birkhoffian associated to such a circuit is not 
regular and not even conservative for nonlinear inductors and capacitors. 




Figure 2 



We have k = 3, p = 2, S/ = 2, Sy = 1, n = 4, m = 4 , b = 8. We choose the reference node to be V5 and 
the current directions as indicated in Figure 2. We cover the associated graph with the loops I\, I2, h, I a- 
Let V = (Vi, V2, Va, Va) G R 4 be the vector of node voltage values, 1 = (i a , i a , Is, , is v ) G R 3 xR 2 xR 2 xR 1 
be the vector of branch current values and v — (v a , v a , vs t , Sy) G R 3 x R 2 x R 2 x R 1 be the vector of 
branch voltage values. 

The branches in Figure 2 are labelled as follows: the first, the second, and the third branch are the 
inductor branches Li, L2, L3, the forth and the fifth branch are the capacitor branches Ci, C2, the next 
two branches are the current source branches S j , Sj 2 , and the last branch is the voltage source branch 



'V, 



The incidence and loop matrices, B £ 9Jts4(R) and A G 9Jls4(R), write as 



B = 



t 


-f 


1 








\ 




-1 
























-1 












-1 












-f 












f 








-1 






f 













\ 





-f 


f 





/ 



A = 



( 
1 

-1 



f 



V 



1 

-1 



-1 1 





\ 



1 

-1 




1 / 



For linear inductors and capacitors, the governing equations have the form: 

f -h-l2+l Sll (t) + l Sl2 (t) = 

-13 - l SIl (t) = 

^-^+v SIl =0 
~^-v sl2 =0 

+ % +V* I2 =0 

§r - §f + (*) = 

where c a 0, a = 1, 2 and h a ^ 0, a = 1, 2, 3, are distinct constants. 



(5.30) 



(5.31) 



Note that i s , , I s 



sv 1 



are given functions of time which describe the currents associated to the in- 



dependent current sources Sj 1 , Si 2 and the voltage associated to the independent voltage source Syi , 
respectively. 
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Once we know the unknowns ii, 12, 13, Qi, Q2, we can determine all the other circuit variables. 
From the first set of equations l|5.31|l . we have 



SVl " dt 



and from the second set of equations l|5.31[l . we conclude 

dl 2 



-L2 



dt 



Therefore, the system 114. If) has now the form 



-II - 19 + 1 SI (t) + h l2 (t) = 

_ dQ 2 _ rfQx _ n 
11 dt dt ~ u 

-I 3 - i sil (i) = 



Ll 



dt C 2 



with 



-1-10 
1 0-1-1 
0-10 



1-10 1 
1-1 



The relations (|3.11() . 14.5(1 . read as follows for this example 

dQ( a ) 



l„ := 



dt 



, a = 1,2,3 



x 1 := Q (1) , x 2 := Q (2 ), x 3 := Q (3) , ,x 4 



1 1 


-1 


1 



Lj x '■— 0,2 



(5.32) 



(5.33) 



(5.34) 



(5.35) 
(5.36) 

(5.37) 
(5.38) 



Using the first set of equations l|5.34|l . we define the 2-dimensional affine-linear configuration space M c 
We solve the corresponding equations l|4.4|l in terms of 2 variables. In view of the notations (|5.37|l . 
we obtain, for example, 



x 2 = -x 1 + f 2 (t) + const 

x 3 = f 3 (t) + const 

r 14 

x = x — x + const 



(5.39) 



with f 2 (t) = J (i Sl (t) + l S2 (t)) dt, f 3 (t) = — J i Sl (t)dt and the other components of / in l|4.6|l being zero. 
Thus a coordinate system on M c is given by 



q 1 := x 1 , q 2 :— x 4 . 



(5.40) 



and the matrices of constants 9t : 



m< 3 

m? 






I 


= 1,2, a = 
3 = 1,2 


( 


1 


\ 




-1 


















1 


\ 


1 


- 1 J 



and C are 



C = 



1 
1 



(5.41) 



23 



In terms of the coordinates 15.40fl . we define the Birkhoffian uj tc = Qidq 1 + Q2dq 2 , as in l|4.8l) - (|4.11ll . that 
is, 



Li +L 2 )q 1 -La- ' 



c 2 c 2 



const 



c 2 



Q 2 (t,q,q,q) = h 1 <? + tW (*) + consi 



Cl c 2 



(5.42) 



Because there exists a loop which contains only capacitors and independent voltage sources, the Birkhof- 



fian (|5.42(1 is not regular. Indeed, the second row of the matrix 



dQ 3 

dq' 



therefore, det 



9Qj 

W 



W=l,2 



contains only zeros, 



0. 



V =1 -. 2 

Though there exists in the network a cutset formed by inductors and independent current sources, 
the Birkhoffian (|5.42|) is conservative in the sense of definition 1(2.19|1 . The function E^it, q, q) is given 
by (j4.13|l . that is, 



-EU {t, q, q) 



-L2 



d 2 f 2 (t) 
dt 2 



+ consti^j q 1 + (v SVl (t) + const 2 ) q 2 



(5.43) 



In order to obtain a regular Birkhoffian we could use the second equation from (|5.42(l and reduce the 
configuration space M c to a vector space M c of dimension 1. The procedure is the same as in the first 
example with linear devices. 

For nonlinear inductors and capacitors, in the coordinate system (15.40(1 on the configuration space 
M c of dimension 2, the Birkhoffian w tc = Qidq 1 + Q 2 dq 2 , where 



Ql(t,q,q,q) 



Q2(t,q,q,q) 



Lx{q l ) 



df 2 {t) 
dt 



df 2 (t)\ d 2 f 2 (t) 



dt 



dt 2 



C 2 (q l — q 2 + const) 

Ci{q 2 ) - C 2 (q 1 -q 2 + const) + v s „ (t) 



(5.44) 



The Birkhoffian (|5.44l) is not regular and not conservative. Indeed, two of the necessary conditions 

^2 ^2 dE OE 

for the existence of the function E UJt : TM — > R such that J2j=i Qj^pi 9> <7i q)dq 3 = Y^j=i dqS ^ 8$ ^ ' 
are 



dq 1 



= + d W) + C,tf - q 2 + const) 



(5.45) 



We easily see that for almost all values of the parameters, d g ^t ^ ^ J^ f 



d 9 V 



0. 



Let us now consider a network that has the oriented connected graph as in Figure 2 in which we 
interchanged the inductor branch L3 with the capacitor branch Ci and the inductor branch L 2 with the 
capacitor branch C 2 . We will see that the Birkhoffian associated to this network is conservative even 
if the inductors and the capacitors in the network are nonlinear devices. 
The system (14.111 has now the form 



ii - 12 - 13 = 

ii(ii)^ + £ 2 (i 2 )^-C 2 ( Q2 ) = o 

-L 2 (i 2 )^ + L 3 (i 3 )^+v Sv (t) = 



(5.46) 
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and 

i SVl = is (5.47) 

v Sll = Ci(q x ) - C 2 (Q 2 ) 

w., a = -C 2 (Q 2 ) (5.48) 

Using the same procedure as above we get the configuration space A4 C of dimension 2. In view of the 
notations l|5.37[l . (|5.38|) . a coordinate system on Ai c is given by 

q 1 := x x ,q 2 := x 3 . (5.49) 

and the Birkhoffian cJt c = Qidq 1 + Q 2 dq 2 , where 

Ql(t, g, <j, g) - [Li^ 1 ) + L 2 ~ q 1 ~ L 2 (q 1 - q 2 ) q 2 - C 2 {-q l + f(t) + const) 
Q 2 (t,q,q,q) = -L 2 (q 1 - q 2 ) q 1 + [L 2 (q 1 - q 2 ) + L 3 (q 2 )q 2 ] q 2 + v SVi (t) (5.50) 

with f b (t) = j (l Sl (t) + i S2 (t)) dt and the other components of / in 1)4.6(1 being zero. The Birkhoffian 
(|5.5UI) is conservative in the sense of definition l|2.19|l . with the function E Ut (t, q, q) given by 



E^faq,*) = I L 1 (q 1 )q 1 dq 1 + J L^q 1 - q 2 ){dq X - dq 2 ) + I L 3 (q 2 )q 2 dq 2 + 
l' 2 (q)(q 1 -<f)dq 1 dq 2 + J J L^dqUq 2 - 
C 2 (-q l + f{t) + cons^dq 1 + V SVl (t)q 2 (5.51) 



where L' 2 :=^0^-. U 
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